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PREFACE 


This  Memorandum  presents  algorithms  forming  the  mathematical 
basis  for  lines-of-communication  (LOC)  models  developed  by  the  author 
and  Captain  M.  J.  Ondrasek,  AFXOXS  (Director  of  Operations,  Systems 
Division,  Hq  USAF) ,  while  members  of  a  special  RAND-AFGOA  (Operations 
Analysis,  Hq  USAF)  joint  Southeast  Asia  study  group. 

The  Memorandum  is  a  generalization  of  earlier  work  by  the  author 
[6];  it  differs  in  that  damage  to  an  arc  resulting  from  interdiction 
may  be  partial  instead  of  total. 

A  model  based  on  the  material  of  this  Memorandum  is  currently 
being  applied  by  RAND  and  AFGOA;  results  will  be  reported  to  the  Air 
Force  as  part  of  a  follow-on  effort  to  the  joint  group's  activities. 
This  model  and  the  results  of  its  application  wil_  be  the  subject  of 
further  Memoranda. 
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SUMMARY 


This  Memorandum  presents  two  general  algorithms  for  targeting  strikes 
in  a  lines-of-communication  (LOC)  network.  Each  depicts  a  situation  in 
which  the  TOCs  are  represented  by  a  network  of  nodes  and  directed  arcs. 
Nodes  are  points  where  two  or  more  arcs  intersect.  Arcs,  which  can  rep¬ 
resent  road,  rail,  or  waterway  segments,  or  points  of  transshipment  be¬ 
tween  two  modes  of  transport,  are  characterized  by  a  beginning  and  ending 
node,  upper  and  lower  bounds  on  flow,  and  a  cost  that  is  a  function  of 
flow. 

The  situations  depicted  by  the  algorithms  differ  only  in  the  form 
of  the  cost  function  on  arcs.  In  the  first  case,  arc  costs  are  linear 
functions  of  flow;  in  the  second  case,  arc  costs  are  convex-piecewise- 
linear  functions  of  flow  with  one  break-point  (i.e.,  two  pieces). 

It  is  assumed  that  the  user  of  the  LOCs  is  attempting  to  achieve  a 
circulation  flow  at  minimum  cost.  This  is  a  very  general  goal  that  in¬ 
cludes,  as  special  cases,  maximizing  flow  between  two  points,  meeting 
a  required  flow  between  two  points  at  minimum  cost,  and  combinations  of 
these  two.  The  algorithms  presented  attempt  to  make  such  costs  as  large 
as  possible  over  time. 

The  effect  of  targeting  strikes  is  to  increase  the  arc  cost  func¬ 
tions  and  to  decrease  arc  capacities  for  a  given  period  of  time.  Thus 
for  the  first  case,  where  arc  costs  are  linear,  the  slope  of  an  arc's 
cost  function  increases  and  its  capacity  decreases  as  the  number  of  strikes 
Increase.  In  the  second  case,  an  arc's  capacity  and  break-point  (i.e., 
the  point  where  the  slope  of  the  cost  function  changes)  are  both  decreas¬ 
ing  functions  of  the  number  of  strikes,  while  the  cost  slopes  before 
and  after  the  break-point  both  increase  with  the  number  of  strikes. 

When  a  single  strike  is  to  be  targeted  against  the  network,  the  arc 
selected  is  one  of  maximum  strike  value,  where  an  arc's  strike  value  is 
defined  as  the  cost  to  repair  it  plus  the  product  of  the  time  to  repair 
it  and  the  cost  increase  of  a  minimum-cost  circulation  flow  resulting 
from  it.  The  strike  value  of  a  particular  arc  may  easily  be  found  by 
finding  a  minimum-cost  circulation  flow  and  its  cost,  first  when  all 
arcs  are  assigned  unstruck  costs  and  capacities,  and  second  when  struck 
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values  are  assigned  to  arcs  under  consideration,  and  then  performing 
the  indicated  arithmetic.  Thus  an  arc  of  maximum  strike  value  could 
be  found  by  calculating  the  strike  value  of  all  arcs.  The  algorithms 
of  this  Memorandum  avoid  this  enumeration  by  taking  advantage  of  the 
fact  that,  when  calculating  the  strike  value  of  one  arc,  upper  bounds 
for  the  strike  values  of  other  arcs  can  often  be  obtained  that  are  low 
enough  to  preclude  them  from  further  consideration. 

In  both  cases,  when  multiple  strikes  are  targeted  against  the  LOCs, 
optimal  arc  selection  is  approximated  by  repeated  application  of  the 
one-strike  algorithm. 
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I .  INTRODUCTION 


This  Memorandum  presents  two  algorithms  for  targeting  strikes 
against  the  arcs  or  links  of  a  combat  force's  lines  of  communication 
(LOCs).  Locations  for  the  strikes  are  chosen  on  the  basis  of  the  degree 
to  which  they  hinder  the  combat  force  or  LOC  user  in  achieving  his  goal. 
This  goal  has  several  possible  formulations.  Among  these  are  maximizing 
flow  between  two  points  and  meeting  a  given  flow  at  minimum  cost. 

The  LOCs  are  represented  by  a  network  of  nodes  and  directed  arcs. 
Arcs  can  represent  road,  rail,  or  waterway  segments,  or  points  of  trans¬ 
shipment  between  two  modes  of  transport.  Each  is  characterized  by  a 
beginning  node,  an  ending  node,  upper  and  lower  bounds  on  flow,  and  a 
cost  that  is  a  function  of  flow.  For  the  situation  treated  by  the  first 
algorithm,  the  cost  function  is  linear;  for  that  treated  by  the  second 
algorithm,  cost  is  convex-piecewise-linear  with  one  break-point.  Thus 
cost  functions  can  be  defined  by  unit-flow  costs.  In  the  first  case 
the  unit-flow  cost  is  constant.  In  the  second  case  it  is  double-valued; 
that  is,  each  arc  has  one  unit  cost  for  flow  up  to  a  certain  value,  and 
another  higher  cost  for  flow  units  beyond  that  value.  These  flow  costs 
could  represent  dollars,  vehicle-hours,  manpower  units,  or  any  other 
appropriate  measure.  Among  quantities  flow  could  represent  are  vehicles 
per  hour  or  tons  per  hour. 

Capacities  (upper  bounds  on  flow)  and  cost  functions  on  arcs  are 
assumed  to  depend  on  the  number  of  strikes  directed  against  them..  For 
the  most  part,  capacities  decrease  and  costs  increase  with  the  number  of 
strikes.  Specifically,  for  the  first  case,  capacities  are  nonincreasing 
and  unit-flow  costs  are  nondecreasing  functions  of  the  number  of  strikes. 
For  the  second  case  both  total  capacity  and  the  point  at  which  unit-flow 
cost  changes  are  nonincreasing  functions  of  strikes,  while  both  unit 
cost  functions  are  nondecre.asing  in  the  number  of  strikes. 

In  addition  to  flow  bounds  and  costs,  arcs  have  repair  times  and 
repair  costs  that  are  associated  with  targeting.  Each  of  these  is  a 
nondecreasing  function  of  the  number  of  strikes.  Thus  the  effect  of 
striking  an  arc  is  to  reduce  its  capacity  and  increase  its  unit-flow 
cost  (or  costs)  for  a  period  of  time  during  which  the  user's  cost  of 
operation  may  increase  and,  in  addition,  to  cause  him  to  incur  a  repair 


cost. 
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It  is  assumed  that  the  general  goal  of  the  user  of  the  LOCs  is  to 
achieve  a  minimum-cost  circulation  flow.  Section  II  will  show  that  this 
includes,  as  special  cases,  maximizing  flow  from  one  node  to  another, 
meeting  a  required  flow  between  two  nodes  at  minimum  cost,  or  combina¬ 
tions  of  these.  The  algorithms  presented  here  attempt  to  allocate 
strikes  against  the  LOC  arcs  in  such  a  way  that,  over  time,  this  minimum 
cost  is  as  great  as  possible.  Such  allocation  is  done  on  the  basis  cf 
immediate  user  cost,  repair  times,  and  repair  cost.  Specifically,  if  a 
single  strike  is  to  be  targeted,  it  is  directed  against  an  arc  of  maxi¬ 
mum  strike  value,  where  strike  value  is  defined  as  the  repair  cost  plus 
the  resulting  cost  increase  of  a  minimum-cost  circulation  flow  multiplied 
by  the  repair  time.  While  this  arc  could  be  selected  by  calculating 
the  strike  value  of  all  arcs,  this  burden  is  av  by  taking  advantage 

of  the  fact  that  in  calculating  the  stri'  r  one  arc,  we  can 

sometimes  obtain  upper  bounds  on  the  stri  _s  of  other  arcs.  We 

can  then  drop  from  consideration  any  arcs  for  which  these  upper  bounds 
are  less  than  the  calculated  strike  values  of  other  arcs  since  they 
could  never  be  arcs  of  maximum  strike  value.  For  multiple  strikes,  no 
method  for  allocating  strikes  optimally,  other  than  complete  enumeration, 
is  known.  It  is  approximated  here,  however,  by  repeated  application  of 
the  one-strike  algorithm. 
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II.  CIRCULATION  FLOWS  ANT)  USER  GOALS 


Finding  a  minimum-cost  circulation  flow  in  a  network  is  done  effi¬ 
ciently  by  Fulkerson's  "Out  of  Kilter"  algorithm  (3).  The  formulation 
of  such  a  problem  is  specified  by  a  network  consisting  of  points  called 
nodes  and  of  directed  arcs,  each  of  which  joins  two  nodes.  An  arc.  orig¬ 
inating  at  node  i  and  terminating  at  node  j  is  designated  by  the  symbol 
(i,j).  Each  arc  (i,j)  has  an  upper  bound,  u^ ,  and  lower  bound,  l^, 
on  itr.  flow,  and  a  cost,  j  •  representing  its  cost  per  unit  of  flow. 
Letting  x^  represent  the  flow  on  arc  (i,J),  the  problem  may  be  formu¬ 
lated  as  follov’a.  Find  flows,  ».  0,  that  minimize 

(1)  I  cij-ii 

ij 


subject  to 

(2)  ^ 

O)  I-jt-I-tj  -°  •“  *• 

)  i 

Expression  (1)  represents  the  total  cost  of  the  flow  pattern,  (2) 
states  that  the  flow  on  each  arc  must  be  between  its  lower  and  upper 
bounds,  and  (3)  rtates  that  the  flow  into  any  node  must  equal  the  flow 
out  of  it. 

Note  that  an  undirected  arc  between  1  and  j  with  sero  lower  bound 
may  be  represented  by  two  directed  arcs,  (l,j)  and  (J,l),  and  that  if 
its  cost  is  positive,  no  minimum-cost  flow  pattern  will  have  strictly 
positive  flow  on  both  of  them. 

Several  possible  goals  for  a  user  of  an  LOC  network  can  be  formu¬ 
lated  as  minimum-cost  circulation  flows  by  proper  data  manipulation. 
The  remainder  of  this  section  discusses  sosm  of  than. 


MAXIMIZING  FLOW 


This  problem  requires  the  maximization  of  flow  from  a  node  that 
shall  be  called  the  source,  to  one  called  the  sink,  subject  to  capacity 
or  upper-bound  constraints  on  the  arcs.  Denoting  the  source  by  S  and 
the  sink  by  T,  the  problem  becomes  one  of  finding  x.  0  that  maximize 
the  source-sink  flow,  v,  subject  to 


(4) 

(5) 

(6) 

(7) 


xij  uij  a11  ^ 


K  ■  \  XS1  + 


v  -  0 


l 

i  j  j 


■  ° an  *  *  s-t 


Zx.T-  l 

i  ]  i 


IT  -  «T,  -  V  '  °- 


If  one  were  to  add  an  artificial  arc  (T,S)  to  the  network  defined 
by  (4)  through  (7),  then  any  aource-alnk  flow  on  the  original  network 
could  be  represented  by  a  circulation  flow  in  the  modified  network. 
Specifically,  (T,S)  would  be  assigned  a  flow  of  v,  and  all  other  arcs 
flows  Identical  to  those  in  the  original  network.  Maximizing  flow  would 
merely  be  maximizing  x^g.  Thus  the  flow  maximization  problem  can  be 
converted  to  a  minimum-cost  circulation  problem  by  constructing  an 
artificial  arc  (T,S)  with  c^s  ■  -1  and  uTg  ■  M,  where  M  is  sufficiently 
large.  All  other  arcs  have  zero  cost,  and  all  arcs  have  ■  0. 

MINIMIZING  COST  OF  REQUIRED  FLOW 


This  problem  specifies  that  a  required  source-sink  flow  be  met  at 
minimum  cost  subject  to  capacity  constraints  on  the  arcs.  If  the  con¬ 
stant  r  is  the  required  flow,  then  this  becomes  a  problem  of  finding 
x^j  x.  0  that  minimizes 


a 

M 


Is  sufficient  since  (4)  and  (5)  imply  v  < 


l  u 

1.3 


U' 
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I  Cij  Xij 
ij 


Xij  ^  Uij  a11  i,;1 

I  xjs  "I  xsj  “  -r 
J  J 

1  Xji  "  I  Xij  “  0  311  1  *  S’T 
j  j 

I  XjT  -  1  *Tj  "  r 

j  J 

This  can  be  converted  to  a  minimum-cost  circulation  flow  problem 
by  adding  an  artificial  arc  (T,S)  with  tTg  »  uTg  ■  r  and  cTg  »  0,  and, 
of  course,  setting  all  other  -  0. 

COMBINING  REQUIRED  AND  MAXIMUM  FLOW  PROBLEMS 

Several  problems  that  are  really  combinations  of  the  maximum-flow 
and  mlnlmum-cost-of -required-flow  problems  may  be  formulated  as  minimum- 
cost  clrc  latlon  flows.  This  section  describes  one  in  which  the  user 
actually  does  desire  to  meet  a  required  flow,  r,  at  minimum  cost.  It 
is  not  known,  however,  whether  the  network  will  support  a  flow  of  r; 
hence  a  more  desirable  formulation  would  be  as  follows.  Find  a  minimum- 
cost  flow  pattern  for  a  required  source-sink  flov,  r,  if  such  a  flow  is 
possible.  If  a  flow  of  r  cannot  be  generated,  find  a  maximum  flow. 

This  can  be  expressed  as:  find  v,  x^  >.  0,  that  minimises 

(13)  ^  ClJ*iJ  ”  ^V* 

IJ 


(8) 

subject  to 

(9) 

(10) 

(11) 

(12) 
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where  M  is  sufficiently  large,  subject  to  v  <  r  and  Expressions  (4) 
through  (7). 

This  converts  to  a  minimum-cost  circulation  flow  problem  by  assign¬ 
ing  the  artificial  arc  (T,S)  with  u^g  =  r  and  c^g  ■  -M,  and  then  setting 
all  lower  bounds  on  arc  flow  to  zero.  After  an  appropriate  pattern  is 
found,  one  may  wish  tc  subtract  out  the  cost  incurred  by  the  flow  on  the 
artificial  arc,  since  this  cost  is  merely  an  artifice  to  Insure  that 
either  a  maximum  flow  or  a  flow  of  r  is  obtained. 

Note  that  if  a  flow  of  r  cannot  be  obtained,  and,  in  addition,  there 
are  several  that  maximize  flow,  this  formulation  will  choose  a  least- 
cost  one  among  these. 


M  ■  Z  c.  .u  ,  is  sufficiently  large. 

l.J  1J 


1 
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III.  ALGORITHMS  FOR  LINEAR  COSTS 

When  arc  costs  are  strictly  linear  functions  of  flow,  strikes  are 
allocated  against  them  on  the  basis  of  upper  and  lower  flow  bounds, 
unit-flow  costs,  repair  times,  and  repair  costs.  All  of  these  except 
the  lower  flow  bounds  vary  with  the  number  of  strikes.  Symbolically, 
these  Inputs  are  written  as  follows: 


(i,j)  designates  the  arc  originating  at  i  and 
terminating  at  j . 


s 


lower  bound  on  arc  flow  for  (i,j). 

upper  bound  on  arc  flow  for  arc  (i,j)  if  s 
strikes  are  directed  against  it. 

cost  per  unit-flow  for  arc  (i,j)  if  s  strikes 
are  directed  against  it. 

repair  time  for  arc  (i,j)  if  s  strikes  are 
directed  against  it. 

repair  cost  for  arc  (i,J)  if  s  strikes  are 
directed  against  it. 


Parameters  with  zero  superscripts  are  essentially  "unstruck" 

parameters  and  normally  all  t?.  and  r*?  would  be  equal  to  zero.  The 

8  ^  8  8  8 
u^j  would  be  nonincreasing  functions  and  the  ,  t^  ,  and  r^  nondecreas¬ 
ing  functions  of  s.  Otherwise,  there  are  no  restrictions  on  these  quan¬ 
tities. 


There  are  several  possible  relationships  for  the  manner  in  which 
these  parameters  vary  with  the  number  of  strikes.  One  can  be  derived 
by  assuming  all  strikes  allocated  against  a  particular  arc  (i,j)  are 
directed  against  the  same  target  and  each  has  an  Identical  and  indepen¬ 
dent  probability  P^  of  successfully  destroying  it.  For  this  situation, 
real  arc  upper  flow  bounds,  unit-flow  costs,  repair  times,  and  repair 
costs  could  take  on  one  of  two  values,  depending  on  whether  the  arc 
was  interdicted  or,  equivalently,  the  target  successfully  destroyed. 


s  s  s  s 

The  u^ ,  c^ ,  ,  and  r^  would  all  be  expected  values  depending  on  the 

probability  of  interdiction.  Letting  the  superscripts  INT  and  UNINT 
refer  to  "interdicted"  and  uninterdicted"  values,  these  expected  values 
would  be: 


s 


u 

c 

t 

r 


ij 

s 

ij 

s 

ij 

s 

ij 


a  -  Cl  -  PtJ 
a  -  a  -  Ptj 

<1  -  (1  -  Py 

(1  -  (1  -  Py 


,a,  INI 

)  )Uy 

Ns.  INT 

)  )Cy 

S  INI 
)  ;ci:} 

.s.  INT 

)  )ry 


+ 
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(1  -  Py 

(1  '  Plj 


. s  UNINT 

>  uu 

. s  UNINT 

>  c« 


For  all  relationships,  however,  allocation  of  strikes  is  done  by 
repeated  use  of  a  one-strike  selection  algorithm.  The  one-strike 
algorithm  selects  for  targeting  the  arc  that  maximizes  repair  cost  plus 
the  product  of  the  repair  time  and  the  cost  increase  of  a  minimum-cost 
circulation  flow.  The  remainder  of  this  section  is  devoted  to  the 
development  of  the  one-strike  algorithm,  the  extension  to  its  use  for 
multiple  strikes,  and  matching  the  inputs  to  conform  to  the  user  goals 
discussed  in  Sec.  II. 


ONE-STRIKE  ALGORITHM 


When  only  one  strike  is  to  be  allocated  against  the  LOCs,  arcs  can, 
of  course,  be  struck  only  zero  or  one  time.  Thus,  for  simplification 
and  for  other  reasons  that  will  become  clear  in  the  next  section,  super¬ 
scripts  for  parameters  dependent  on  the  number  of  strikes  will  be  elimi¬ 
nated  when  referring  to  the  one-strike  algorithm  and  a  prime  will  be 
used  to  denote  struck  values.  These  are  now  written  as  follows: 


'ij 

ij 

ij 

ij 


ij 


upper  bound  on  arc  flow  for  (i,j)  if  unstruck, 
upper  bound  on  arc  flow  for  (i,J)  if  struck, 
cost  per  unit-flow  for  arc  (i,j)  if  unstruck, 
cost  per  unit-flow  for  arc  (i,j)  if  struck, 
repair  time  for  arc  (l,j). 
repair  cost  for  arc  (i,Jn 
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For  arc  selection,  let  k  be  the  cost  of  a  minimum-cost  circulation 
flow  if  no  arcs  are  struck,  and  k^  the  cost  of  a  minimum-cost  circula¬ 
tion  flow  if  arc  (i,j)  is  struck  but  all  other  arcs  are  unstruck.  Then 
the  value  of  striking  arc  (i,j)  is  defined  as 


(14) 


v(i, j )  *  (k^  -  k)t  +  r  , 


The  one-strike  algorithm  finds  the  arc  of  maximum  strike  value.  This 
arc  could  be  found,  of  course,  by  actually  calculating  the  strike  value 
for  each  arc.  However,  the  procedure  for  calculating  the  strike  value 
of  an  arc  (or,  more  specifically,  the  calculation  of  k^)  yields  a  cir¬ 
culation  flow  that  can  sometimes  be  used  to  calculate  upper  bounds  on 
strike  values  of  other  arcs.  Of  course,  an  arc  with  such  an  upper  bound 
that  is  less  than  the  calculated  strike  value  of  some  other  arc  cannot 
be  an  arc  of  maximum  strike  value. 

More  specifically,  suppose  the  calculation  of  a  minimum-cost  cir- 

» 

culation  flow  when  arc  (a,b)  is  struck  yields  an  arc  flow  x^j  —  uij  ^or 

some  ( i  ,j ) .  Such  a  flow  would  be  feasible  when  arc  (i,j)  is  struck,  and 

totaling  up  the  arc  costs  for  this  flow  pattern  would  yield  an  upper  bound 

for  k^j  .  This  upper-bound  calculation  would  differ  from  the  calculation 

of  k  only  in  that  the  struck  cost  is  used  for  arc  (i,1)  instead  of 
ab 

for  arc  (a,b),  which  is  now  assigned  its  unstruck  cost.  Thus  one  obtains: 


(15) 


k  <  k  .  + 
ij  -  ab 


(c 


ij 


-  C  )x  .  -  (c  -  c  .  )x  , 
ij  ij  ab  ab  ab 


Substituting  the  right-hand  side  of  Expression  (15)  fcr  k^  in  (14) 
yields : 


(16)  v(i,J)  <  (kab  +  (ClJ  -  ClJ)xn 


c  »)*  u 

ab  ab 


■t 


Expression  (16)  will  be  used  ir  the  one-strike  algorithm.  Before 
the  algorithm  is  presented,  several  terras  need  to  be  defined  formally, 
some  of  which  have  been  discussed: 
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k  ■  cost  of  minimum-cost  circulation  flow  if  no  arcs 
are  struck. 


cost  of  a  minimum-cost  circulation  flow  if  arc 
(i,j)  is  struck  and  all  others  are  unstruck. 


v(i,J)  -  the  value  of  striking  arc  (i,j)  as  defined  by  (14). 


upper  bound  on  the  value  of  striking  arc  (i,j). 


w  -  the  arc  of  highest  strike  value  among  all  those 
whose  strike  valuas  have  been  calculated. 


Note  that  if  v(w)  2.  v^  for  all  (i,j),  then  w  is  the  arc  to  strike. 

The  algorithm  will  start  with  w  undefined,  v(w)  ■  0,  and  all  v^  ■  ». 

Then  it  will  find  a  minimum-cost  circulation  flow  with  no  arcs  struck, 

and  lower  the  v  in  accordance  with  a  slight  modification  of  Expression 
ij  i 

(15)  wherever  i.  u^  Then  k^  and  v(a,b)  are  calculated  for  a 

particular  arc,  (a,b),  which  has  the  property  that  v  ,  >.  v  for  all 

ab  lj 

(i,J).  The  quantity  v#fc  la  then  lowered  to  v(a,b)  and,  when  possible, 
other  v^  values  are  lowered  in  accordance  with  the  right-hand  side  of 
Expression  (16),  thus  bringing  them  closer  to  the  v(l,j)  values.  If 
v(a,b)  2.v(w),  then  (a,b)  replaces  w.  Otherwise,  w  remains  unchanged. 
Then  (a,b)  Is  replaced  by  another  arc  of  highest  v^  and  the  procedure 
is  repeated  until  a  point  is  reached  where  v(w)  2.  v^  for  all  (i,J). 

At  this  point,  w  is  the  arc  of  maximum  strike  value  and  v(w)  its  strike 
value.  Formally,  the  algorithm  is  as  follows. 


1.  Set  v(w)  »  0  and  all  v^  -  •. 

2.  Using  unstruck  costs  and  upper  bounds  for  arcs,  find  a  minimum 
cost  circulation  flow  and  denote  its  value  by  k.  For  all  arcs  with  x 


iJ 


i  u^ ,  set 


iJ 


(c 


IJ 


”  CiJ)xiJtlJ  *  riJ* 
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3.  Choose  an  arc,  (a,b),  such  that  v^  .  If  =  v(w), 

terminate,  since  w  is  the  arc  to  strike  and  its  strike  value  is  v(w). 
Otherwise  continue. 

t 

4.  Change  (a,b)'s  upper  flow  bound  to  u  ,  and  its  unit-flow  cost 

I 

to  c  .  and  find  a  minimum-cost  circulation  flow.  Denote  its  value  by 
ab 

kab‘ 

5.  Set  vab  -  (kab  -  k)tab  +  rab.  If  vab  v(w),  set  w  -  (e,b) 

and  v(w)  *  v  .  .  Otherwise,  leave  w  and  v(w)  unchanged. 

lb  t  • 

6.  For  each  arc  with  x^  <.  u^  ,  calculate  A^  M  (k^  +  -  c^) 

x  -  (c'  -  c  )  x  -  k}t..  +  r.  and  if  A  <  v  ,  set  v  =  A... 


'ab  ab  ab 


■ij  iJ 


ij  ij 


ij  ij 


(This  step  may  be  eliminated  for  arcs  with  v^  i.v(w),  since  such  arcs 
have  essentially  been  eliminated  as  candidates  to  strike). 

7.  Restore  the  upper  flow  bound  of  (a,b)  to  u^  and  the  unit  flow 
cost  of(a,b)  to  c  ^  and  return  to  step  3. 

Step  1  is  merely  an  initialization  step.  Step  2  calculates  k  and 
attempts  to  lower  the  v .  Step  3  tests  to  see  if  the  arc  of  maximum 
strike  value  has  been  found  (i.e.,  if  v(w)  v^  for  all  (i,j))  and  if 
not,  selects  another  arc  for  examination.  Steps  4  and  5  find  the  value 
of  striking  the  arc  being  examined.  Step  5  also  tests  to  se  if  the  arc 
being  examined  has  a  strike  value  higher  than  those  of  all  other  arcs 
that  have  previously  been  examined  and,  if  so,  designates  it  as  such. 
Step  6  attempts  to  lower  the  v^  of  all  arcs  with  v^  >  v(w),  and  step  7 
restores  the  parameters  for  the  examined  arc  to  unstruck  values. 


MULTIPLE  STRIKES 


There  are  no  known  methods,  other  than  complete  enumeration,  for 
allocating  strikes  optimally,  when  several  strikes  are  to  be  targeted. 

Thus  optimality  will  be  approximated  by  allocating  strikes  one  at  a  time. 
Each  time,  the  arc  against  which  to  allocate  the  next  strike  is  deter¬ 
mined  by  the  one-strike  algorithm.  In  applying  the  one-strike  algorithm, 
arc  parameters  are  determined  by  the  number  of  strikes  already  allocated 
to  them.  Spec i f ica 1 lv ,  if  the  first  n  strikes  include  s  to  be  directed 

against  arc  (1,1),  then,  in  applying  the  one-strike  algorithm  to  allocate 
—  St 

the  n  +  1  strike,  the  unstruck  parameters  for  arc  (i,J)  correspond  more 
or  less  to  those  for  s  strikes  and  the  struck  parameters  to  s  +  !  strikes. 
More  specifically,  they  are 
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uij 

m 

Uu 

1 

s+l 

uiJ 

m 

uij 

s 

cij 

cij 

1 

s+l 

ciJ 

m 

cij 

S+1 

ij 

cij 

s+l  s 

rU 

m 

riJ  -rij 

The  generalized  algorithm  for  allocating  n  strikes,  then,  is  as 
follows : 


t 


ij 


1. ^  For  each  archer  -  uj.,,  -  ujr  cy  -  cir  , 

■  t^j,  and  r^j  »  r  -  r^ .  Set  fi  -  1. 

2.  Use  the  one-strike  algorithm  to  find  the  location  for  the  n1*1 


strike. 


3.  Let  (a,b)  be  the  location  for  the  nth  strike  and  let  s  be  such 
8 

that  u  .  -  u  ,  (i.e.,  among  the  first  n  -  1  strikes,  s  is  the  number 

ab  ab  « 


allocated  to  (a,b)).  Then  set  u 


8+1  s+2 

^  'ab  "  uab  ’  ab  uab  ’ 

cs+2  .  ts+2  .  r»+2  _  rs+l 

ab  *  ab  ab  ab  ab  ab  * 


'ab 


s+1 

'ab 


'ab 


4.  If  n  -  n,  terminate.  All  strikes  have  been  allocated  and  the 
number  of  strikes  allocated  to  (i,j)  is  s(l,j),  where  s(i,j)  satisfies 
u^  •  •  Otherwise,  Increase  n  by  one  and  return  to  Step  2. 


Step  1  initializes  parameters.  Step  2  finds  the  location  for  the 
nth  strike.  Step  3  updates  parameters  for  the  arc  selected.  Step  4 
tests  to  see  if  all  strikes  have  been  allocated. 


USER  GOALS  AND  INPUTS 


Section  II  shewed  that  several  possible  user  goals,  while  not  appear¬ 
ing  to  consist  of  finding  a  minimum-cost  circulation  flow,  could  be 
formulated  as  such.  Three  such  goals,  each  of  which  Involved  finding 
a  source-sink  flow,  required  the  construction  of  an  artificial  arc  from 
the  sink  (T)  to  the  source  (S).  User  goals  were  determined  by  the  values 
of  the  upper  and  lower  flow  bounds  and  unit-flow  cost  for  (T,S).  Since 
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sorae  of  these  parameters  depend  on  the  number  of  strikes,  their  values 

must  be  generalized.  Furthermore,  val.ies  must  also  be  assigned  to 

repair  times  and  repair  costs  for  (T,S). 

The  generalization  of  the  parameters  for  the  artificial  arc  are 

such  that  its  strike  value  is  always  zero.  Thus  repair  costs  and  repair 

times  are  zero  for  all  levels  of  strike  and  equal  to  those  defined  in 

s  s 

Sec.  II.  Specifically,  in  all  three  special  cases,  t^,,  -  0  and  rTg  ■  0 
for  all  s.  If  the  LOC  user  is  attempting  to  maximize  flow  from  S  to  T, 

HTS  -  0,  u®s  -  M,  and  c®g  -  -1  for  all  s,  where  M  is  very  large.  If  the 

g 

goal  is  to  meet  a  required  flow,  r,  at  minimum  cost,  Sl^g  ■  r,  u^g  ■  r, 

3 

and  c^,g  *  0  for  all  s.  In  the  combined  problem  where  the  user  is  attempt¬ 
ing  to  meet  a  required  flow,  r,  at  minimum  cost  if  such  a  flow  can  be 
achieved,  but  to  maximize  flow  if  a  flow  of  r  cannot  be  achieved,  £„  ■ 

S3  J  ^ 

0,  uTg  ■  r,  and  cTg  ■  -M  for  all  s,  where  M  is  again  very  large. 
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IV.  ALGORITHMS  FOR  PIECEWISE  LINEAR  COSTS 


* 

A  RAND-AFGOA  study  group  found  it  desirable,  where  arc  parameters 
are  based  on  expected  values  that  in  turn  are  based  on  probabilities 
of  successful  targeting,  to  let  unit-flow  costs  reflect  the  fact  that 
arcs  which  are  little  used  when  unstruck  may  not  be  used  at  all  when 
successfully  targeted.  While,  ideally,  one  might  wish  to  allow  for 
arc  costs  to  be  strictly  convex  functions  of  flow,  the  desirability 
of  using  the  programmed  "out-of-kilter"  algorithm  led  to  piecewise- 
linear  cost  functions  with  one  break-point. 

This  requires  new  definitions  for  the  upper  flow  bounds  and  unit- 
flow  costs  for  arcs.  These  are: 


upper  bound  on  arc  flow  for  arc  (i,j)  if  s  strikes 
are  targeted  against  it. 


—  8 

u^  value  at  which  unit-flow  cost  changes  on  arc  (i,j)  If 
s  strikes  are  targeted  against  It. 

g  _  g 

c^  unit-flow  cost  for  units  up  to  u^  on  arc  (i,j)  If  s 

strikes  are  targeted  against  It. 

8  —8 
c^j  unit-flow  cost  for  units  in  excess  of  u^  on  arc  (i,j) 

If  8  strikes  are  targeted  against  it. 


''.he  parameters  ,  t®j ,  and  r*j  are  as  defined  when  arc  costs  are 
linear  functions  of  flow.  The  cost  of  a  circulation  flow  pattern  with 
flows  would  (with  superscripts  deleted)  be 


(17) 


l  ^ 


X  + 

ijij 


i 


C  0  + 

ijij 


I 


:lj  (XU 


V* 


XU>UiJ 


It  Is  assumed  that  u 


‘ij’  1.1  ’ 


and  r 


s 

ij 


8  .  S  —  g 

and  G  are  nonincreasing  in  s,  and  c^ , 
are  all  nondecreasing  in  s.  Of  course,  u^  <.  u®j 


See  Preface. 


-15- 


and  cij  cij  • 

In  Che  RAND-AFGOA  study  mentioned  earlier,  it  is  assumed  that  all 
strikes  directed  against  arc  (i,J)  are  aimed  at  the  same  target  and  each 
has  independent  probability  p^  of  interdiction  or  successful  targeting. 
The  upper  flow  bound  and  unit-flow  cost  on  arc  (l.j)  are,  respectively, 
u^  and  c^  when  interdicted,  and  u^  and  c^  when  unlnterdlcted . 
Arc  parameters  as  a  function  of  the  number  of  strikes  are  expected  values 
obtained  when  it  is  assumed  that  arcs  not  used  to  full  capacity  when  un¬ 
interdicted  will  not  be  used  at  all  when  Interdicted.  Thus,  they  are 


u 


ij 

_s 

UU 


'IJ 

s 

ij 


,,  .  s  UNINT  .  ...  INT 

•a-Vuu  *  °  '  °  '  "u*  ’“ij 

.s  UNINT 

■ (1  -  V  "ij 

UNINT 

m  r 

ij 

_  INT 

U 

■  (l  -  (l  -  p.j)  )ttJ 


a  -  a 


INT 


pu,‘)ru 


INT  INT 

where,  of  course,  t^  and  r^  are  repair  times  and  repair  costs  when 
arc  (1,J)  is  interdicted. 

As  in  the  case  of  linear  costs,  strikes  are  allocated  one  at  a  time 
by  a  one-strike  algorithm  that  chooses,  as  the  arc  to  strike,  the  one 
that  maximizes  repair  cost  plus  the  product  of  repair  time  and  the  cost- 
increase  of  a  minimum-cost  circulation  flow  pattern. 


*  —  f  8 

Since  c^  <_  c^j,  the  problem  of  finding  a  minimum-cost  circulation 

flow  can  be  converted  to  one  in  which  arc  costs  are  linear  functions  of 
flow.  This  is  done  by  replsclng  each  arc,  (1,)),  by  two  arcs  In  parallel. 
The  first  of  these  has  upper  flow  bound  of  u*^  and  unit-flow  cost  of  c*^ , 

while  the  second  has  values  of  u|^  -  u*^  and  for  these  parameters. 

These  inputs  are  now  suitable  for  the  "out-of-kllter"  algorithm. 
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ONE-STRIKE  ALGORITHM 


As  in  the  esse  of  linear  costs,  an  Iteration  of  the  one-strike 
algorithm  will  consist  of  finding  the  strike  value  of  some  arc,  and 
then  trying  to  lnprove  upper  bounds  already  obtained  for  the  mike 
value  of  the  arcs  whose  strike  values  have  not  been  calculated.  The 
algorithm  terminates  when  the  calculated  strike  value  of  some  arc  is 
at  least  as  great  as  the  upper  bound  on  strike  value  for  any  other  arc. 
The  notatlonal  simplification  adopted  earlier  of  dropping  superscripts 
and  designating  struck  values  by  a  prlae  will  also  be  adopted  here. 
They  are  now  written: 


C!J 

*U 

rU 


upper  bound  on  arc  flow  for  (l,j)  if  unstruck 

upper  bound  on  arc  flow  for  (i,j)  if  struck 

value  at  which  unit-flow  cost  changes  on  arc  (i.J)  If  unstruck 

value  at  which  unit-flow  cost  changes  on  arc  (i.J)  if  struck 

unit-flow  cost  for  units  up  to  u^  on  (i.J)  if  unstruck 

unit-flow  cost  for  units  up  to  on  (i.J)  if  struck 

unit-flow  cost  for  units  in  excess  of  u^  on  (i.J)  if  unstruck 

unit-flow  cost  for  units  in  excess  of  u^  on  (i.J)  if  struck 

repair  time  for  (i.J) 
repair  cost  for  (I.J) 


The  quantities  k,  k^  ,  v(i,J),  v^ ,  and  w  are  as  defines  In  the  one- 
strike  slgorlthai  for  linear  costs. 

If  a  ■lnlaum-cost  circulation  flow  obtained  when  (a.b)  Is  struck 
yields  x^  ltijj,  then  this  flow  is  feasible  when  (1,))  is  struck  and 


i  k  .  ♦ 

ab 


1J 


-  5 


ab 


(18) 


where 
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61J  "  (clj 


c1j)(«in[u1j .*tj )) 


♦  (c^  -  c^ ) (ux[0,  alnfx^.u  )  -  Ujjl) 

♦  (cjj  -  c  HaaxfO.x^  -  8  ]). 

Th«  quantity  6  la  the  lncraaaa  In  coat  on  arc  (1,J)  obtained  by  sub¬ 
stituting  struck  paraaeters  for  unatruck  ones.  The  first  tens  In  the 
right-hand  side  of  (19)  represents  the  cost-increase  on  (i,J)  for  flow 

_•  i 

units  up  to  u£j»  the  second  for  flow  units  between  0^  and  u^,  and  the 
third  for  flow  units  in  excess  of  u^  .  Substituting  the  right-hand  side 
of  (18)  for  k^  In  (14),  the  expression  for  v(l,j),  one  obtains 


v(l.J)  i  (k#b  ♦  -  6^b  -  k)!^  ♦  rtj . 


Thus,  the  one-strike  algorlthai  for  finding  the  arc  of  aaxlnun  strike 
value  Is  as  follows: 

1.  Set  v(w)  •  0  and  all  •  •. 

2.  Using  unstruck  costs  and  upper  bounds  for  all  arcs,  find  a 

alnlaua-cost  circulation  flow  and  denote  Its  value  by  k.  For  all  arcs 
with  Xjj  4.  u^,  act  v^  •  where  4^  ia  as  defined  in  (19). 

BAX 

3.  Choose  an  arc,  (a,b),  such  that  v^b  -  ^  j)vij‘  **  Vab  "  v^* 
terninate,  since  w  is  the  arc  to  strike  and  its  strike  value  Is  v(w). 
Otherwise  continue. 

4.  Change  (a,b)'s  upper  flow  bound  to  u^,  the  point  at  which 

unit-flow  cost  changes  to  0  b»  and  Its  low  and  high  unit-flow  costs  to 

£^b  and  c^b  respectively.  Then  find  a  ninlaiuai-cost  circulation  flow 

and  denote  Its  value  by  k  . 

ab. 

5.  Set  v#b  •  (k-b  -  k)t-b  ♦  r#b.  If  v^  2.  v(w)  set  w  -  (a.b) 
and  v(w)  •  v^b>  Otherwise  leave  w  and  v(w)  unchanged. 

6.  For  each  arc  with  x  t.  u  ,  calculate  A  •  (k  .  ♦  4  - 

l  J  X  j  X  j  AD  X j 

«ab  ■  k>*n  4  rij*  *nd  lf  *‘tj  *  vij*  ,et  vij  ■  Aij‘  <Thl*  st*P  “*y  be 
eliminated  for  arcs  with  v  1  v(w).) 
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7.  Restore  (a,b)'s  upper  flow  bound  to  u^*  the  point  at  which 
unit-flow  cost  changes  to  u  .  ,and  low  and  high  unit-flow  costs  to  c  . 

u  D  3D 

and  c^  respectively,  and  return  to  step  3. 

The  functions  performed  by  the  individual  algorithm  steps  are 
essentially  the  same  as  those  in  the  linear-cost  case.  Specifically, 
step  1  initializes;  step  2  calculates  k  and  attempts  to  lower  the  v^ ; 
step  3  tests  to  see  if  the  arc  of  maximum  strike  value  has  been  found 
and,  if  not,  selects  another  one  for  examination;  steps  4  and  5  find 
the  value  of  striking  the  arc  being  examined  and  if  this  value  exceeds 
v(v),  lets  it  replace  w;  step  6  attempts  to  lower  the  v^ ;  and  step 
7  restores  the  parameters  of  the  examined  arc  to  unstruck  values. 

MULTIPLE  STRIKES 


The  optimal  allocation  of  several  strikes  against  the  LOCs  is  ap¬ 
proximated  by  repeated  application  of  the  one-strike  algorithm,  with  arc 
parameters  updated  for  each  application.  If  the  one-strike  algorithm  is 
entered  with  s  strikes  already  targeted  to  arc  (i,J),  then  the  input  pa¬ 
rameters  for  arc  (i.j)  are 


8+1 

‘ij 
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ij 
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Following  the  foraac  for  Che  linear-cost  case,  Che  general  algo¬ 
rithm  for  targeting  n  scrikes  in  the  pleeewise-llnear-cost  case  is  as 


follows : 


0  '  1 

1.  For  each  arc,  sec 

-0  -'  -1  O’  1  „ 

Cij*  Cij  "  Cij  *  Cij  "  Cij’  CiJ  "  CiJ*  CiJ 


uir  au  ■  a°r  5i)  ■  v 


hr  *nd  ru 


■  r 


1J 


See  H  •  1. 


2.  Use  Che  one-scrike  algorithm  to  find  the  location  for  Che 

-th  ... 
n  strike. 


3.  Let  (a,b)  be  the  location  for  Che  strike  and  let  s  be 

s  s+l  *  s+2  -  _  -s+l  „ ' 

such  that  uab  -  u4b.  Then  set  uab  -  u-b  ,  u#b  -  «*ab  •  u 

s+l  '  s+2 

Cw"cv.C.»C.  , 

ab  ab  ab  ab 


ab 


s+2  _ 

u  .  ,  c  . 

ab  ab 


-  c8+1 


«'  _  rS+2 
ab  ’  ab  ab  1 


t  -  t8*1 

u  u 


ab  1  ab 
and 


s+2 

r  ,  •  r  . 
ab  ab 


-  r 


8+1 

ab 


4.  If  n  -  n,  terminate.  All  strikes  have  been  targeted  and 
the  nuaber  of  strikes  targeted  to  (i,J)  is  a(i,J),  where  s(i,J)  satisfies 
Otherwise,  increase  fi  by  one  and  return  to  step  2. 


USER  GOALS  AND  INPUTS 


With  respect  to  the  special  user  soals  discussed  in  Sec.  II, 
parameters  for  the  artificial  arc,  (T,S),  are  quite  similar  to  those 
of  the  linear-cost  case.  All  u8^  ■  0  and  C®^  may  take  on  arbitrary 
values.  This  essentially  makes  cost  on  the  artificial  '.re  a  strictly 
linear  function  of  flow,  with  c*<,  the  unit-flow  cost  v  :n  s  strikes 
have  been  targeted  against  it.  Thus  the  newly  defined  parameters  have 
essentially  been  zeroed  out  for  this  arc  and,  for  all  user  goals  dis¬ 
cussed,  the  quantities  r8^  ,  t8^,  u” j •  *nd  take  on  the  sane 

values  as  in  the  linear  cost  case.  Specifically,  t*g  ■  0  and  r*g  •  0 
for  all  s.  When  the  LOC  user  is  attempting  to  maximize  flow  from  S  to 
T,  t^g  •  0,  u^!,.  ■  M,  and  c8g  •  -1  for  all  s;  when  attempting  to  meet  a 
required  flow,  r,  at  minimum  cost,  l^g  ■  r,  u8g  •  r,  and  c8g  •  0  for 
all  s;  and  when  attempting  to  meet  a  required  flow  r,  at  minimum  cost 
if  possible,  and  If  not  possible,  to  maximize  flow,  t^g  -  0,  u*g  ■  r, 
and  c8g  •  -1  for  all  s.  In  all  of  these  cases  M  represents  a  very  large 
number. 
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